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We obtain exact travelling wave solutions for three families of stochas- 
tic one-dimensional nonequilibrium lattice models with open boundaries. 
These solutions describe the diffusive motion and microscopic structure 
of (i) of shocks in the partially asymmetric exclusion process with open 
boundaries, (ii) of a lattice Fisher wave in a reaction-diffusion system, 
and (iii) of a domain wall in non-equilibrium Glauber-Kawasaki dynam- 
ics with magnetization current. For each of these systems we define a 
microscopic shock position and calculate the exact hopping rates of the 
travelling wave in terms of the transition rates of the microscopic model. 
In the steady state a reversal of the bias of the travelling wave marks a 
first-order non-equilibrium phase transition, analogous to the Zel'dovich 
theory of kinetics of first-order transitions. The stationary distributions 
of the exclusion process with n shocks can be described in terms of n- 
dimensional representations of matrix product states. 

PACS numbers: 05.70.Ln, 82.20.Mj, 82.40.Fp, 02.50. Ga 
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1 Introduction 



Systems of diffusing and reacting particles are usually described on the 
macroscopic level by hydrodynamic equations for coarse-grained quantities 
like the particle density which represent the order parameter specifying the 
macroscopic state of the system ^Ij. Paradigmatic examples for these equa- 
tions are the Burgers equation for driven diffusive systems with particle con- 
servation 12] or the Fisher equation for reactive systems without conservation 
law 13 El- These equations are in general non- linear and exhibit shocks in 
some cases. This means that the solution of the macroscopic equations may 
develop a discontinuity even if the initial particle density is smooth. In or- 
der to understand the emergence of such behaviour from the microscopic 
laws that govern the stochastic motion and interaction of particles it is nec- 
essary to derive the macroscopic equations from the microscopic dynamics 
rather than postulating them on phenomenological grounds. To solve this 
problem it is evident that detailed insight in the microscopic structure of 
non-equilibrium systems exhibiting macroscopic discontinuities must be ob- 
tained. 

A considerable body of results of this nature has been obtained for specific 
one-dimensional lattice models defined on the integer lattice 2 El HI , the 
best-studied example being the asymmetric simple exclusion process (ASEP) 
[HI El • In this basic model for a driven diffusive system each site k is either 
empty (n^ = 0) or occupied by at most one particle (n^ = 1). A particle on 
site k hops randomly to the site k + 1 with rate Dr and to the site k — 1 
with rate Di, but only if the target site is empty. Otherwise the attempted 
move is rejected. The jumps occur independently in continuous time with 
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an exponential waiting time distribution. In the hydro dynamic hmit the 
system is described by the Burgers equation which exhibits shocks. Such 
a shock discontinuity may be viewed as the interface between stationary 
domains of different densities. Relaxing the requirement of particle number 
conservation leads to class of systems which are generically called reaction- 
diffusion processes, but which by a variety of mappings also serve as toy 
models for non-conservative spin-flip dynamics (in the context of magnetic 
systems), epidemic spreading, growth processes and transport phenomena in 
biological and ecological systems and elsewhere HO] • 

Most of the results for the dynamical behaviour have been obtained for 
infinite particles systems. In many of the physical applications, however, 
one has to study finite systems with open boundaries where particles are 
injected and extracted. This is crucial to take into account as - in the 
absence of equilibrium conditions - the boundary conditions determine the 
bulk behavior of driven systems, even to the extent that boundary induced 
phase transitions between bulk states of different densities occur [HI [T21 CB] ■ 
Qualitatively, the strong effect of boundary conditions on the bulk can be 
attributed to the presence of steady-state currents which carry boundary 
effects into the bulk of the system. Quantitatively, exact results for the 
steady state of the ASEP have helped to show that part of the nonequilibrium 
phase diagram of driven diffusive systems with open boundaries, viz. phase 
transitions of first order, can be understood from the diffusive motion of 
shocks ^lE]; in analogy to the Zel'dovich theory of equilibrium kinetics 
of first-order transitions. As in equilibrium, the nonequilibrium theory of 
boundary-induced phase transitions requires the existence of shocks which 
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are microscopically sharp. 

In a series of recent papers JHI - jlHI these considerations, originally for- 
mulated for conservative dynamics, have been extended to non-conservative 
react ion- diffusion systems. Moreover, there are exact results about shocks 
in reaction-diffusion systems with branching and coalescence j2n]-|2S] (here 
shocks are known as Fisher waves on the macroscopic scale) and in spin- 
flip systems where shocks correspond to domain walls 26J. However, no 
exact results have been reported so far for non- stationary travelling waves 
in open systems. Here we wish (i) to establish a complete picture about 
exact travelling wave solutions for the specific family of systems to which 
these processes belong (viz. single-species exclusion processes with two-body 
nearest neighbour interaction and no internal degrees of freedom) and (ii) to 
study the dynamics and microscopic structure of these travelling shocks in 
systems with open boundaries. Since many of the powerful techniques used 
for treating the ASEP do not apply to nonconservative systems we propose 
a general approach that can be applied to any lattice model: We take as ini- 
tial distribution a shock distribution with given microscopic properties and 
determine the class of models for which the shock distribution evolves into a 
linear combination of similar distributions with different shock positions. In 
this paper we identify three families of processes with this property. 

The paper is organized as follows: In the following section we define the 
class of models that we consider and we also define shock measures for these 
systems. In Sec. 3 we determine the families of models with travelling wave 
solutions on the finite lattice. This is followed by some new results for the 
ASEP with open boundaries in Sec. 4. In Sec. 5 we summarize our results 
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and draw some conclusions. 

2 Reaction-diffusion systems and shock mea- 
sures 

2.1 Stochastic single-species models 

We consider Markovian interacting particle systems of a single species of 
particles without internal degrees of freedom which have hard-core two-body 
interactions with their nearest neighbour sites. We describe hard-core inter- 
action due to excluded volume in terms of an exclusion processes where each 
lattice site may be occupied by at most one particle. This class of models 
may therefore be described by a set of occupation numbers n = {rii, . . . , til} 
where = 0, 1 is the number of particles on site A; on a lattice of L sites. 
There is a one-to-one correspondence to classical spin systems where the oc- 
cupation number rik = represents spin up while Uk = I represents spin 
down. 

The stochastic dynamics are defined in terms of transition rates (transi- 
tion probabilities per infinitesimal time unit). The process is fully defined 
by the 12 rates Wij for changes of the configuration of a pair of neighbouring 
sites k and k + 1 j2ZI : 



where i = 1,2,3,4 is the decimal value plus one of the target configuration 
(^fc'^fc+i) read as a two-digit binary number and j is the respective value of 
the initial configuration (n^n/j+i), as shown below. 




(1) 



Diffusion to the left and right (01 ^ 10) W32, W23 
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Coalescence to the left and right (11 10,01) W34, W24 

Branching to the left and right (10,01 11) W43, W42 

Death to the left and right (10, 01 00) W13, W12 

Birth to the left and right (00 — * 10,01) W31, W21 

Pair Annihilation and Creation (11 ^ 00) Wu, w^i 

From time to time we also use the more intuitive symbols Dr = W23, Di = W32 
for the hopping rates. Notice that combinations of individual processes may 
describe other physically meaningful processes. E.g. coalescence and death 
with equal rates is equivalent to single-site radioactive decay (1 — 0) with 
that rate. The inverse of the rate is the mean life time of a particle. For 
injection and extraction of particles at the boundaries we introduce the rates 

Injection and Extraction at the left boundary (0 ^1) ct, 1 
Injection and Extraction at the right boundary (0 ^1) S, f3 

The time evolution is defined by a continuous-time master equation for 
the distribution P{ni, ■ ■ ■ ,nL;t) which we write in terms of the quantum 
Hamiltonian formalism The distribution is mapped to a probability vec- 
tor I P(t) ) which contains as components the probabilities P{ni, ■ ■ ■ , ul', t). 
The time evolution is generated by the stochastic Hamiltonian H whose ma- 
trix elements are the transition rates between configurations. The Markovian 
time evolution can then conveniently be cast in the form of an imaginary time 
Schrodinger equation 

j^\P{t))-H\P{t)) (2) 



6 



with the formal solution 



Pit)) 



-Ht\ 



P{0)). 



(3) 



The quantum Hamiltonian H for the family of processes defined above 
has the structure 

H=J2hk + h + bL. (4) 

Here 



k=l 



hh = - 



( ■ Wyi Wn Wl4 \ 

W2\ ■ W23 W2A 

W31 W32 . W2.A 

\ W41 W42 W43 . / 



(5) 



k,k+l 



is the local transition matrix acting nontrivially on sites /c, /c + 1. The diag- 
onal elements are the negative sum of the transitions rates in the respective 
column, as required by conservation of probability. The boundary matrices 



hi 



-a 7 
a —7 



-5 /3 
6 -13 



(6) 



generate the boundary processes. 

The invariant measures | P* ) of the process, i.e., the stationary probabil- 
ity distributions, satisfy the eigenvalue equation 



H\P*) = 0. 



(7) 



We stress that the analogy to quantum mechanics is a formal one, for details 
see 0. 

The equations of motion for the expected local particle density take the 
form [2Z] 
d 



dt 



+D,{n,^i{t)n,{t) ) + D2{n,{t)n,+i{t) ) 



with 

^1 = W21 + W41 Bi = W23 + W43 - ""^21 - W41 

Cl = W12 + Ws2 + W21 + W41 Di = Ci - W23 - W43 - Wi4 - W34 

M = W31 + W41 B2 = W32 + Uli2 — W31 — Wii 

C2 = Wis + W23 + W3I + W41 D2 = C2- W32 - W42 - - W24- 

(9) 

In analyzing these equations the question arises how to treat the nonhnear- 
ity in the lattice equation, i.e., the two-point correlator Di{nx^i{t)nj.{t) ) + 
D2{nx{t)nx^i{t)) . Calculating its time- derivative leads to a coupling to 
three-point correlation functions and eventually to a hierarchy of equations 
which is just as untractable than the master equation itself. Only for some 
families of models the system of equations decouples and exact results can be 
obtained 1771 |2Hl 12^1 ISO] • Therefore here we do not follow this traditional 
approach but rather investigate the time evolution of the measure. 

2.2 Product measures and shock measures 

The stationary distribution of the family of processes defined above depends 
on all the transition rates and is not known in general. On some parameter 
manifolds, however, the stationary distributions are simple Bernoulli product 
measures 

P*(2l) = n((l-P)'^n.0 + P5n,,l) (10) 

where the probability of finding a particle at each site A; is p and independent 
of the occupation of other lattice sites, i.e., where there are no correlations. 
It is easy to see that P* depends only on the total number = X^fc of 
particles in the configuration n, one has P*{n) = (1 — p)^~^ . 
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In the quantum Hamiltonian formulation this distribution is represented 

by a tensor state 

\P')=[ /) ^\P) (11) 
The family of processes for which this is an invariant measure can be deter- 



mined easily from ((Zj). One first determines the manifold of bulk rates 



1-3 



such that 

hk\ p) = A{nk+i - nk)\ p) (12) 

with an arbitrary constant A. This yields three equations for the 12 bulk 
parameters. The solutions define a manifold of processes with uncorrelated 
stationary distributions, provided the system has periodic boundary condi- 
tions. In order to satisfy ((Tj) for systems with open boundaries one determines 
the boundary parameters by the relations 

hi\ p) = Ani\ p) , hL\p) = -AhlIp). (13) 

For each boundary this is one equation for two rates. Notice that these 
relations contain the stationary particle density p as free parameter. 

Bernoulli shock measures are product measures with a jump in the density 
at some site m (jT]). They are represented by a tensor state 

T) r) • 

There are no correlations, but the density in the left domain of sites 1, . . . , /c 
is pi and then jumps to p2 in the right domain + 1, . . . , L of the system. 
Since there are no correlations one may regard the lattice unit as the intrin- 
sic shock width. Hence shocks which are described by a such a distribution 
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are microscopically sharp and have a very simple internal structure, charac- 
terized by the absence of any correlation between particle positions. There 
is no process of the form Q for which a shock distribution with a shock 
at some given site k is stationary. However, as shown in next section, linear 
combinations of shock measures may be stationary distributions. Notice that 
a shock at position m = corresponds to a Bernoulli measure with density 
P2 while a shock at position k = L corresponds to a Bernoulli measure with 
density pi. 



Pi 



P2 



Figure 1: Density profile of a Bernoulli shock measure 



3 Shocks as stable collective excitations 

A linear combination of shock measures is a stationary measure if a given 
shock measure evolves into a linear combination of shock measures after time 
t., i.e., where | P{t) )m,p^^p2 = ^^P )pi,P2 

I^WL,P„P2 = E^'K^'^)l^)p.P2- (15) 

k=0 

The physical interpretation of this property is that a shock retains its internal 
structure at all times, only the position of the shock is shifted by a random 
amount. The probability of moving from the initial shock position m to site 
k after time t is the quantity p{m, k; t). Hence we shall refer to measures of 
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the form (jl5|) as diffusive shock measures. The random walk nature of such a 
shock gives rise to the interpretation as a collective single-particle excitation. 
This can be made more precise by the implications of (fTH|l . The evolution 
into a linear combination of shock measures implies for an infinitesimal step 
(which is generated by H) the evolution equation 

-H\m = di\m-l)^^^^^+d2\m+l ) ^^^^^ - {di + ds) | m )^^^^^ (16) 

which is the evolution equation for a biased single-particle random walk with 
hopping rate d2 to the right and hopping rate di to the left. Hence we 
have to determine processes such that (fT^ is satisfied. The family of shock 
distributions defined by the densities pi, p2 forms an invariant sector lA under 
the time evolution of the system. Notice that ()16p implies the existence of at 
least two stationary product solutions in a periodic system. These stationary 
states have densities pi,p2 respectively. The boundary rates have to chosen 
such that at the left boundary pi is stationary, while at the right boundary 
P2 is stationary. 

Solving (fTBj) leads to three classes of reaction-diffusion models with an 
invariant sector U which are described below. 

3.1 Asymmetric Simple Exclusion Process (ASEP) 

The simplest process which satisfies (|16p is the ASEP with hopping to the 
left and right (without loss of generality we assume a bias to the right) and 
injection and extraction at both boundaries. Hence the nonvanishing rates 

W32, W23, 7> «, 5. (17) 
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These rates together with the densities pi and p2 satisfy the following con- 
ditions: 

a(l - pi) - 7pi = (w23 

13 p2 - 5(1 - P2) = (^23 

Both densities have to fuUfil the conditions < pi < p2 < 1 in this case. 
Condition p8|) was obtained for the infinite system in [3T]) using symmetry 
properties of the quantum Hamiltonian (see next section). In the bulk of 
the lattice the shock position moves like a biased lattice random walk with 
hopping rates 

d, = {D,.-Di f'^^~P'^ (21) 

P2 - Pi 

to left {i = 1) and right (z = 2) respectively. This leads to the well-known 
exact expressions for the shock velocity Vg = d2 — di and shock diffusion 
coefficient Ds = (^2 + di)/2 O E2] as long as the shock is far away from the 
boundaries. 

The new result are conditions p9|) . ()20|) which imply that at the bound- 
aries the shock is reflected. According to the properties of a biased random 
walk on a finite lattice with reflecting boundaries and bulk rates fl21|) its 
stationary position after equilibration is geometrically distributed, i.e., the 
probability of finding the shock at site k on the lattice is of the form 

p*(A:)cx(^|y. (22) 

Depending on the bias of the shock the steady state of the system is in 
the low-density subphase Aj (for ^2 > di), in the high-density subphase Bj 
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Pi)Pi 



W32 



Pi(l - P2) 
W32)pi(l - Pi) 

«^32)P2(1 - P2) 



(18) 
(19) 
(20) 



(for di > ^2) or on the first-order coexistence line (for di = d2) f2lES]- From 
(|22p we read off the exact inverse localization length 

r' = ln(rf2)-ln(rfi) (23) 



which was conjectured in ^2] to describe the localization of the shock through- 
out the subphases Aj and Bj. In fact, since the existence of a shock is a 
generic property of driven diffusive systems, our results support the picture 
developed in |T3] where it is argued that a localization length of the form 
is universal for driven diffusive systems in the subphases Aj and Bi. 



3.2 Branching— Coalescing Random Walk (BCRW) 

In this case we have the following non-vanishing rates: 

W34, W24, W42, W43, W32, W23, a, 7, P (24) 

In the periodic system there are two translation invariant stationary states 
in this model: Bernoulli measures with zero density and with a density p* 
resp. which depends only on the ratio 

W24 + ^^34 ^ 1 -P* ^25) 

Wa2 + W43 P* 

between the branching and coalescence rates. For the existence of diffusive 
shock measures in the open system one of the densities pi has to be zero, 
the other has to be p* . Without loss of generality we set p2 = 0. The 
non- vanishing rates then have to satisfy the conditions 

W23 = / (26) 

p* 

7 = - — + (1 ^ p*)wz2 - - — r~^43 + p*wm- (27) 

p* p* 
13 



This leaves 7 free parameters. 

Instead of the BCRW we could have chosen a birth-death-diffusion model 
with rates obtained from (f^ - PH|l by interchanging particles and vacancies. 
The densities behave under this transformation according to pj — > 1 — pj for 
i = 1, 2 As microscopic shock position it is convenient to choose the position 
of the rightmost particle. 

From |TH|l we obtain the shock hopping rates 

di = w^s/p* (28) 
d2 = w^sil-p*)/p*. (29) 

Hence the velocity of the shock 

Vs = d2 - di = W43 (30) 

is determined by the rate of branching to the right. The diffusion coefficient 
reads 

D = wj-^ = ^ fl + 2^^^Hi±^) . (31) 
2p* 2 V Wi2 + WiJ ^ ' 

This quantity remains finite in the limit p* — 0. 

For a special tuning of the coalescence rates (instantaneous on-site coales- 
cence which is equivalent to D^. = W24,, Di = W34) this process can be solved 
exactly with the help of the so called inter-particle distribution functions 
(IPDF) [2HIEIIESI5 or, equivalently, using free fermion techniques, reviewed 
in detail in By passing to the continuum limit (lattice constant a ^ 0) 
ben-Avraham has shown for the free-fermion choice of coalescence rates with 
infinitesimal branching rate (proportional to the lattice constant) that the 
model has shock-like solutions if the initial state has zero density on one side 
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and stationary density on the other side of the origin [231 • He has has also 
derived properties of higher order correlations which suggest the existence 
of a diffusive shock measure at least for this special limiting case of rates. 
Thus one has a picture comparable to the situation in the asymmetric sim- 
ple exclusion process, however, with the difference that in the asymmetric 
exclusion process the densities on both sides of the shock are arbitrary while 
in the coalescence-branching model these densities are fixed to be zero and 
p* = pi respectively. Moreover, in the limit of ben-Avraham the stationary 
density p* in the active domain is non-zero, but infinitesimal. Straightfor- 
wardly extending our result to the infinite system proves the existence of 
such a diffusive shock measure and shows that such a shock solution persists 
also for finite densities in the active domain. 

We remark that for the case of symmetric hopping Dr = Di interesting 
macroscopic dynamics arise from (jH)) if we consider infinitesimal rates of 
branching and coalescence of the order and also rescale time by a? (diffusive 
scale). We set 

W24 = a^W24, W34 = a^W34, W42 = a^W42, W43 = a^W43 (32) 
and obtain within mean field theory 

f = P^ + A)0 + Mp*-p) (33) 

with k = W2A + 'W34 + W42 + 1&43. This is the usual Fisher equation [H] which 
is also known to have travelling wave solutions similar to shocks. 

Mean field theory for infinitesimal branching and coalescence rates resp. is 
justified since in this case the dynamics in any finite region are dominated by 
the pure exclusion process and hence is expected to be locally stationary and 
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hence to have no correlations ^Tj. Neglecting correlations in the derivation 
of the hydrodynamic limit of (jH)) for the asymmetric process and keeping 
terms up to second order in the lattice constant one obtains 

with infinitesimal viscosity v = a{Dr + Di), single-particle velocity v = D,.— 
Di + a{wi2 — W43), nonlinearity = 2{Dr — Di) + a{wi2 + — W43 — W34) 
and reaction constant k = W42 + W24 + W43 + ^34- This equation was studied 
by Murray [SE] where it was shown that there are shock-like travelling wave 
solutions. For a = the stationary equation reduces to an ordinary first 
order differential equation. With fixed boundary densities the solution can 
be obtained using the approach of Ref. pTj . 

3.3 Asymmetric Kawasaki-Glauber Process (AKGP) 

In this case the nonzero rates are the death and branching rates as well as 
one hopping rate where without loss of generality we consider nonvanishing 
hopping to the left 

W12, Wis, u!42, W43, ''^32, a, P (35) 

In the absence of diffusion {ws2 = 0) this model is a variant of zero-temperature 
Glauber dynamics [26j with a non- vanishing magnetization current [3T. In- 
cluding diffusion corresponds to a nonequilibrium coupling of the zero-temperature 
process to an infinite-temperature heat bath with asymmetric Kawasaki spin 
exchange dynamics. In a biological context branching corresponds to cell 
duplication by mitosis which can occur only if there is space available for a 
second cell. In this setting the death process describes the effect of certain 
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types of drugs which kill both cells in the event of mitosis [2H] • The two sta- 
tionary densities are and 1 respectively. The diffusive shock measures are of 
the same form as for the ASEP and BCRW resp. with pi = 1, p2 = 0. Given 
an initial step function profile with a single domain wall ...1111100000... it 
is clear that the only events that can occur are the hopping of the domain 
wall to the right or left. Thus the domain wall performs a lattice random 
walk just as in the previous examples with hopping rates 

di = Wis, d2 = W43 (36) 

to left (di) and right (^2) respectively. 

We remark that for a special choice of the branching rates W43 = w^ — Di, 
W42 = Wis + Di the nonlinear contributions to the equations of motion (jH)) 
vanish identically [27^. Hence the exact evolution of the density profile is 
given by a lattice diffusion equation for all initial distributions. In this case 
the equations of motion for the density p{x, t) do not give any indication of 
the existence of stable shocks. Hence the existence of a non-linearity in the 
dynamical equation for the density is not necessary for having shocks in the 
associated process. 

4 ASEP with open boundaries 

Here we wish to explore some of the ramifications of the results of the previous 
section on shock diffusion in the ASEP with open boundaries. 

1) As discussed above the stationary distribution of the shock position de- 
scribes the steady state properties of the exclusion process with open 
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boundaries along the manifold of boundary parameters defined by (|19|) - 
(j2ni)- The exact steady state properties of the exclusion process for 
all values of the boundary parameters may be calculated explicitly ei- 
ther by solving recursion relations ^2] or by using the so-called ma- 
trix product approach [[3] which involves the representation theory 
of a quadratic algebra equivalent to a g-deformed harmonic oscillator 
algebra jHSI- The conditions (fT^ - POjl are equivalent to the condi- 
tions for the existence of a two-dimensional representation of the Fock- 
representation of the quadratic algebra used to calculate the stationary 
state properties of this model in Refs. [TBlEnilini- This can be seen as 
follows: Let us define the function Kj. as: 



-X + y + W23 - W32 + J {-X + y + W23 - W32y + 4:xy 
f^+{x,y) = 

(37) 

Then the conditions ()18 p - ()20|) can be written in the form: 

Pi = T—^f T (38) 

1 + K+(7,a) 

= TWM) '''' 

W23 = n+{j,a)K+{6,P)w32 (40) 
We remind the reader that higher dimensional representations satisfy 

{W23/W32T = i^+h, (y)K+{S, p), (41) 

for the derivation see [10] • 

2) Diffusive shock measurees have also been considered for the infinite 
system. This was done in for discrete time evolution (parallel 
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updating) and in [^I] for continous time evolution. The mathemat- 
ical structure behind the one-particle nature of the bulk shock mo- 
tion is the t/g[S'f/(2)]-symmetry of the generator of the ASEP with 
reflecting boundaries which relates one-particle states to shock states 
through the action of the ladder operator of Ug[SU{2)]. Surprisingly 
the open boundary conditions considered here break the symmetry, 
Yet the single-particle nature of the shock remains. This is reminiscent 
of partially broken symmetries observed in spin chains with diagonal 
boundary fields |12]. 

3) In the infinite system consecutive multiple shocks which each satisfy 
(fTHj) evolve according to n-particle dynamics, i.e., the Uq[SU (2)]-symmetry 
relates these n-shock measures to states with n particles |3T]. Using 
the ansatz discussed above it is easy to verify that by imposing the 
boundary conditions ()4H) required for the existence of n-dimensional 
representations of the stationary quadratic algebra one also obtains 
closed equations of motion for shock measures in the open system. To 
see this one adopts a slightly different definition of the shock mea- 
sure introduced in [21]: It is a product measure with density 1 at the 
shock positions ki and intermediate densities pi between sites ki-i,ki. 
The consecutive densities each satisfy ()18|) which by iteration leads to 
the condition ()4H) for the existence of n-dimensional representations of 
the stationary algebra. Hence the n-dimensional representations of the 
stationary algebra describe the stationary linear combination of shock 
measures with n consecutive shocks. Notice that this modified defi- 
nition of shock measures also allows for a representation of arbitrary 
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shock measures (not satisfying any special relation between consecutive 
densities) in terms of linear combinations of the special shock measures. 

4) The definition of the shock position by the jump property of the shock 
measure (or the presence of a particle with probability 1 at the shock 
position in the alternative definition, respectively) is not a microscopic 
definition of the shock position for a single realization of the process. 
As the system evolves in time one cannot trace the shock position to 
tell where exactly the shock is located. In a single realization the shock 
position would emerge only after spatial coarse graining. A convenient 
definition of the microscopic random position Xt of the shock in a single 
realization of the process is the position of the second class particle 
E1I33|. This particle behaves like an ordinary particle with respect to 
empty sites and like an empty site with respect to an ordinary particle. 
The second class particle has a drift towards the shock j3] and its 
position may thus be defined as the position of the shock. (Cf. [46j for 
this choice.) Its diffusion coefficient has been obtained in Ref. [221 • The 
density profile of the invariant shock measure as seen from the second 
class particle was calculated in jlTj. Its shape depends on the hopping 
ratio q = ^Jd^/Di and on the limiting densities. For limiting densities 
satisfying ()18p the exact state is a pure Bernoulli shock measure with 
densities pi to left of the second-class particle and p2 to its right. This 
observation suggests investigating the dynamics of shock measures with 
second class particles at the shock positions. 

We apply again the strategy of following the time evolution of shock 
measures by calculating the action of the Hamiltonian on the measure. 
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In the infinite system one finds that indeed these measures form a closed 
sector analogous to lA if the condition (fTH|) is satisfied for consecutive 
shock densities. For the study of the open system we also need to define 
the properties of the second class particle at the boundary sites. To 
this end we denote second-class particles by the symbol B and explicitly 
consider reservoir sites 0, L + 1 which may either contain a -B-particle 
with probability 1 or no i?-particle, but an A-particle with probability 
pi,2 respectively. We denote these two possible configurations of the 
boundary reservoir by -Ri and B respectively and represent the left 
boundary processes as follows: 

i?iO — > RiA with rate SipiD^ 

RiA^RiO with rate - pi)A 

RiB —>■ RiA with rate 6ipiDr 

RiB i?iO with rate - pi)A (42) 

SO -> RiB with rate 5iA 

BA RiB with rate 6iDi. 

At the right boundary we define analogously two reservoir states i?2, B 
resp. on site L + 1 and allow for processes with rates 



0R2 - 


^AR2 


with rate 


S2P2D1 




AR2 - 


^0R2 


with rate 


62(1- 


P2)Dr 


BR2 - 


^AR2 


with rate 


S2P2D1 




BR2 - 


-^0i?2 


with rate 


62(1- 


P2)Dr 


OB - 


^BR2 


with rate 


S2D1 





(43) 
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AB BR2 with rate 820^. 

The number of S-particles is conserved in this dynamics. Physically 
this corresponds to the reflection of shocks at the boundaries in the 
open system. The first two transition rules for both boundaries which 
do not involve 5-particles satisfy the injection/absorption rules ()19|). 

Considering the system with several 5-particles one may study the 
time evolution of n consecutive shocks with increasing densities at each 
point of discontinuity. For multiple shock measures with consecutive 
densities satisfying 

A il-p^+l)pi ^ ' 

one finds n-particle dynamics. On the hydrodynamic scale n consec- 
utive shocks obey simple deterministic n-body dynamics: They move 
with constant speed until two shocks meet and then coalesce into a sin- 
gle shock. Thus after some time only one shock (the leftmost, which is 
the fastest) survives [48 . For the special family of boundary densities 
considered above this phenomenon can be studied on the lattice scale. 

5 Conclusions 

We have studied the dynamics of a class of react ion- diffusion processes with 
open boundaries on the lattice scale and established a complete list of mod- 
els where exact travelling-shock solutions exist. For these systems we have 
detailed knowledge about the microscopic structure of the shock. We found 
that there a three families of such models: The ASEP, the BCRW, but 
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on a more general manifold of parameters as considered previously, and a 
Kawasaki-Glauber spin-flip dynamics. In all three cases the time evolution 
of the shock measure is equivalent to that of a random walker on a lattice 
with L + 1 sites with homogeneous hopping rates in the bulk and special 
reflection rates at the boundary. The existence of such processes imply a 
rather remarkable property. Shocks behave like collective single-particle ex- 
citations already on the lattice scale - not only after coarse-graining where 
all the microscopic features of the shock are lost. The reduction of the expo- 
nentially large number of microscopic internal degrees of freedom (2^) to an 
only polynomially large number of macroscopically relevant degrees of free- 
dom (L + 1) is not an uncontrolled and only phenomenologically motivated 
approximation, but an exact result on all scales of observation. 

As is long known from zero-temperature Glauber dynamics a hydrody- 
namic description of an interacting particle system in terms of a PDE for 
the particle density may fail to give any hint on the microscopic structure of 
the macroscopic solution even if the hydrodynamic equation is exact. Our 
results for the AKGP indicates that this property is not a special feature only 
of Glauber dynamics. It remains as an open question under which general 
conditions and in which way the presence of a stable shock in a stochastic 
interacting particle system is reflected in the structure of the hydrodynamic 
limit. It also would be interesting to investigate travelling shocks in systems 
with defects, where in the case of the ASEP exact results are available for 
the steady state |13 EOl EH El • 
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